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Abstract 

We determine numerically the single-particle and the two-particle spectrum of 
the three-state quantum Potts model on a lattice by using the density ma- 
trix renormalization group method, and extract information on the asymptotic 
S-matrix of the quasiparticles. We find that the finite size spectra can be com- 
pletely understood in terms of a simple effective model introduced in a previous 
work, and are consistent with an asymptotic S-matrix of an exchange form. 
The finite size spectra are, on the other hand, inconsistent with a diagonal 
asymptotic S-matrix, predicted by the bootstrap solution of the corresponding 
perturbed conformal field theory. 



1. Introduction 

Being the simplest generalization of the transverse field Ising model, the q 
state quantum Potts model is one of the most paradigmatic models in statisti- 
cal physics and quantum field theory. The case of q = 3 is somewhat peculiar 
and is of particular interest. In one spatial dimension, the q = 3 state quan- 
tum Potts model displays a second order quantum phase transition between a 
ferromagnetic state and a paramagnetic state, just as the transverse field Ising 
model P, 0, Q ■ The properties of the critical state itself are very well charac- 
terized: at the critical point, an exact solution is available and the scaling 
limit is known to be described in conformal field theory (CFT) by the minimal 
model of central charge C = 4/5 [H], Q, with the so-called D4 partition function 
in the ADE classification Q. The ordered and disordered phases of the Potts 
model are, on the other hand, much richer than those of the transverse field 
Ising model: Similar to e.g. antiferromagnetic chains of integer spins j^. [oL Hoj|. 
the gapped phases (i.e., the ferromagnetic as well as the paramagnetic phase) 
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possess excitations with internal quantum numbers; as a consequence, the dy- 
namics of these quasiparticles are much richer than those of the transverse field 
Ising model. 

While the critical properties of the Potts model are very well understood, 
the quasiparticle properties of the gapped phases are somewhat controversial. 
In the continuum limit, the properties of the Potts model are thought to be 
described by the so-called scaling Potts field theory, which is a perturbation of 
the fixed point conformal field theory, uniquely determined by the symmetries. 
The application of the machinery known as the S-matrix bootstr ap II I f yields a 
diagonal quasiparticle S-matrix for low energy particles [l2|, EH, Il4, Il5l. Hil . [l7j 
and implies that the internal quantum numbers of two colliding particles are 
conserved during a scattering process (see fig. [l]a) i.e. the scattering process 
should be pure transmission. The consistency of the bootstrap S-matrix with 
the perturbed conformal field theory is demonstrated in a separate paper by 
one of the authors [3] . 

However, recent calculations on the lattice version of the quantum Potts 
model contradicted these results Q: perturbative calculations in both phases 
as well as rather strong renormalization group arguments yielded a coherent 
picture, and supported that in the gapped phases, rather than being diagonal, 
the asymptotic S-matrix of the lattice Potts model assumes the " universal" form 
emerging in various spin models [§], as well as in the sine-Gordon model fill ] : 
S — > —X, with X the exchange operator. According to this result, quasiparticles 
of small momenta should scatter on each other by exchanging their quantum 
numbers (see fig.[T]b). We must emphasize that the structure of the asymptotic 
S-matrix has important physical consequences: an S-matrix of the exchange 
form yields universal diffusive finite temperature spin-spin correlation functions 
at intermediate times 0, [Til l2oj |. while a diagonal S-matrix would result in 
exponentially damped correlations [ 1 OL 1 2 lL \22\. 

Although the arguments of Refi Q are very robust, the results of Ref. Q 
were met by some skepticism. We therefore decided to study in detail the 
two-particle spectrum of the q = 3 state quantum Potts model on a lattice 
using the powerful numerical method of density matrix renormalization group 
(DMRG). With this method we are able to compute the finite size spectrum very 
accurately, and compare it to the predictions of the simple effective theory of 
Ref. Q, and to those of the bootstrap. As we shall see, the finite size spectra are 



2 



in complete agreement with the theory of Ref. [3j and an asymptotic S-matrix 
of the exchange form, while they obviously disagree with the predictions of the 
bootstrap, and therefore exclude the corresponding diagonal S-matrix. As we 
discuss later, we believe that this discrepancy between the continuum approach 
and the lattice calculations can be traced back to the (incorrect) assumption of 
intcgrability, and the presence of " dangerously irrelevant" operators. 

2. The Potts model and its quasiparticles 

In its lattice version, the Potts model consists of a chain of generalized 
spins having internal quantum states with i labeling the lattice sites and 
/i = l,...,q the possible internal states of the spins. The Hamiltonian of the 
q-statc quantum Potts model is then defined as 

i fi— 1 i 

Here the traceless operators Pj* = \n)a(fi\ — 1/q tend to project the spin at 
site i along the "direction" /x, and thus the first term of cq. (UJ promotes a 
ferromagnetic ground state, with all spins spontaneously polarized in one of the 
directions, In contrast, the second term in cq. ([I]) represents a "transverse 
field", with the traceless operator Pi = |An)ii(An| — l/q trying to align the 
spins along the direction |Ao) = Ia*)/\/<7- The relative strength of these 
two terms is regulated by the dimensionless coupling, g. These terms obviously 
compete with each other, and their competition leads to a phase transition: 
for large values of g one finds a paramagnetic phase with a unique ground 
state, while for small g a ferromagnetic phase appears with q degenerate ground 
states, spontaneously breaking the global S g symmetry. In the q = 3 case, — 
on which we focus here, — the transition occurs at a coupling g = g c = 1, and 
it is of second order: quasiparticles are gapped on both sides of the transition, 
but the quasiparticle gap A vanishes continuously at the transition as A ~ 
J \g-l\^ 6 

The q state Potts model obviously possesses a global § q permutation sym- 
metry. As a consequence, the global cyclic permutation Z\fi)i = \fx+ 1 mod q)i 
leaves the Hamiltonian also invariant, and can be used to classify its eigenstates 
as 

Z\Q) = e mQ \Q), (2) 

with Q an integer and the angle f2 defined as = 2?:/q^ In the particular case 
of q = 3, considered here, Q can take values of Q = and Q = ±. In this case, 
pairwise spin exchanges (e.g., fj, = 1 -<-> 2) also imply that states with quantum 
numbers Q = ± come in degenerate pairs. 

The structure of quasiparticles in the ferromagnetic {g < 1) and in the 
paramagnetic (g > 1) phases can be easily understood in the perturbative limits, 
g <C 1 and j>1. For g > 1 the ground state |0) is unique, and quasiparticles 
consist of local spin flips of §3 charges Q = ±. For g < 1, on the other hand, the 



1 This holds even in the ferromagnetic phase, but there states with spontaneously broken 
symmetries must be mixed. 
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ground state is 3-fold degenerate, |0) — > |0) p , and quasiparticles correspond to 
domain walls between these ground states, [i — > fjf = /j, + 9 mod 3, with 9 = ± 
the quantum number of the domain wall. 

Similar to the Ising model, the Potts model is known to be self-dual. High- 
temperature - low-temperature duality [23[ in the d = 2 classical Potts model 
implies a duality jf}l/j for the quantum Potts model 24j. In the Appendix 
we show that duality holds even on the level of the matrix elements of the 
Hamiltonian, and therefore one can map the spectra in the Q = sectors for 
g and 1/g by simply rescaling the energies with appropriate factors. We thus 
have 

E$= (g)=gE$= (l/g) (3) 

for all eigenstates n with periodic boundary conditions (PBC), as also verified 
later numerically. This duality relation has important consequences, and shall 
allow us to relate various energy- and length scales on the two sides of the 
transition. 



3. Effective theory and two particle S-matrix 

In an infinite system, the elementary excitations of the gapped phases can 
be classified by their momentum, k, and for small momenta their energy can be 
approximated as 

e(fc) = A+^ + ... (4) 
2m 

independently of their internal quantum number. Here m = m(g) is the quasi- 
particle mass, and A = A(g) denotes the quasiparticle gap. 

In the very dilute limit, interactions between quasiparticles can be described 
in terms of just two-body collisions, and correspondingly, by just two-body 
scattering matrices and interactions. Assuming pairwise and short ranged in- 
teractions between the quasiparticles, one thus arrives at the following effective 
Hamiltonian (in first quantized form) 0, [13] , 

Nqp 1 FP- 

11 » ^7^' • X **) + ••■ - (5) 

i— 1 1 i<j 

with Xi and Ui denoting the coordinates and internal quantum numbers of 
the quasiparticles, and N qp their number. The above Hamiltonian acts on 
many-particle wave functions ip{ai}({%i})i which are bosonic (invariant under 
exchanges (xi,ai) <-> (xj,crj)), and correspond to states of the form = 
J2{x i }J2{a l }' t f ; {<yi}({ x i})\{ x i}A (J i})- Tli e dots in eq. © denote higher order 
terms, which are irrelevant in the renormalization group sense, and do not in- 
fluence the asymptotic low-energy properties of the theory. 

The scattering of two quasiparticles on each other can be characterized by the 
two-particle S-matrix, which, in view of the energy and momentum conservation, 
has a simple structure. The two-particle S-matrix, in particular, relates the 
amplitude of an incoming asymptotic wave function ipkia-i ,ka<ri ( x i ^ x ?) ~ 
A™ x cr ^(ki,k2)e llyklXl+k2X2 } with quasiparticle momenta k\ > k^ to that of the 
outgoing wave function, ipk^nM^^i > x %) ~ b °ctT,t 2 ( fc i> k 2 )e l( - klXl+k2X ^ as 

B out = S(h - fca) A in . (6) 
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The structure of the two-body S-matrix is further restricted by §3 symmetry: 



S(k) 



( s 3 {k) \ 

si(fc) s 2 (k) 

s 2 (k) si(fc) 

V s 3 (k) J 



(7) 



In the following, we shall only investigate the scattering of quasiparticles in the 
channels {-) — } and { — h}. In these channels, the eigenvalues of the S-matrix 
read 

s t (k) = e M *> = Sl (k) + s 2 (k) , (8) 
s s (k) = e 2i5 °W = Sl (k)-8 2 (k) , (9) 

where we introduced the "triplet" and "singlet" eigenvalues, st(k) and s s (k), 
and the corresponding phase shifts, 5t(k) and S s (k). As shown in Ref. |3|, 
interactions in the singlet channel are irrelevant for k — > (the wave function 
has a node at x\ = x 2 ), while they are relevant in the triplet channel, unless 
some very special conditions are met by the effective interactions Q . As a result, 
generically one finds St(k — > 0) = — 1 while s s (k — > 0) = 1, as also confirmed by 
direct calculations in the g — > 00 and g — > limits Q. As a consequence, by 
analyticity, the phase shifts must have the following small momentum expansion: 

s t(k) =~\ s S n ( fc ) +a t k+... , 6 s (k) = -a s k + ... . (10) 

Notice that these expressions (together with ss(k 0) = —1) give rise to a 
low-momentum scattering matrix of the form, S w —X. In contrast, perturbed 
conformal field theory yields a diagonal low- momentum S-matrix with st(k — >• 
0) = 1, corresponding to irrelevant interactions even in the triplet channel. 
However, this would require very special interactions, and is not guaranteed by 
S3 symmetry. 

3.1. Two-particle spectra: paramagnetic phase 

The two-particle spectrum of a finite system of size L a = 1 follows from 
the asymptotic form of the S-matrix. In the following, we shall focus exclusively 
on the simplest case of periodic boundary conditions (PBC). 

In the paramagnetic phase, quasiparticles carry a "chirality" label, a = Q = 
±. Therefore, the Q = + sector of the spectrum contains single quasiparticle 
excitations of chirality a — + [described by eq. flU] as well as, e.g., two-particle 
excitations with charges <j\ = a 2 = — . As a consequence, in the Q = ± sectors 
it is numerically hard to separate two-particle states from the single-particle 
states. We therefore focus on the sector Q = 0, where single quasiparticle states 
are absent, and, above the ground state, the spectrum starts directly with two- 
particle eigenstates of quasiparticles with charges 01 = ± and <j 2 = =p. 

For large system sizes, the quantization of the momenta k\ and k 2 of the 
quasiparticles is determined by the periodicity condition on the wave function, 
$f(xi,x 2 ) = ^(xi + L,x 2 ) = ^(xi,x 2 + L), and, just as in Bethe Ansatz, the 
energy is the sum of the two quasiparticle energies, E = e(ki) + e(k 2 ). Taking 
particle i = 1 around the system (see fig. l3.1[) then yields the following condition, 



C\ \ _ ikl L ( si(ki-k 2 ) s 2 (ki-k 2 )\f Ci 
C 2 J e \ sa^-fe) Sl (h-k 2 ) J\C 2 



(11) 
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Figure 2: (Color online) Translating one of the quasiparticles around by system size L in the 
paramagnetic (left figure) and ferromagnetic (right figure) phases. 



with C\ = A+_(ki > k 2 ) and C 2 = A™, (fci > k 2 ) the wave function amplitudes 
for < x\ < x 2 < L, defined earlier. Taking particle i = 2 around, x 2 — > x 2 +L, 
yields a similar equation. In the triplet channel, C± = C 2 , we thus obtain 



st(ki 



-ik\ L 



-ik 2 L 



(12) 



Using the asymptotic expansions of the phase shifts, eq. (JTUJ), and solving 
cq. (fT2"j) to leading order in 1/L then gives 



ni ,n 2 



2 A = £ 



1 , 1 (rii — n 2 + l) 2 , 
-(ni + n 2 2 + -— 7 ' + Oil L ) 



where n\ and n 2 denote integers, and we introduced the energy unit, 

1 /2tt n 



(13) 



(14) 



In cq. (|13[) . to comply with the bosonic nature of the excitations, the quantum 
numbers n\ and n 2 must satisfy ni > n 2 . 

The previous analysis can be carried over to the singlet sector, C\ = —C 2 , 
with little modification, and there it yields the following finite size spectrum: 



Ei 



-2A = £ 



1 , 1 (Til — 7l2) 2 , 9s 

4f"' + " 2 ) + 4 (i-4 +0(1/£) 



(15) 



However, now ri\ and n 2 must satisfy ni > n 2 since for rti = n 2 the wave 
function vanishes trivially. 



3.2. Two particle spectra: ferromagnetic phase 

As discussed earlier, the ground state of the infinite system in the ferromag- 
netic phase has broken §3 symmetry, and correspondingly, it is 3-fold degenerate, 
|0) M . □ Excitations are kinks (domain walls), and the corresponding two-particle 
states read 

\x\ Oi,x 2 2 )i_ l , (16) 



2 Here we used brackets rather than angular brackets, to explicitly emphasize that the states 
10)^ are interacting many-body cigenstates of the Hamiltonian. 
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with \i the vacuum polarization at x — > — oo, X\ the positions of the domain 
steps, and 6i = ± the step sizes. As we shall also demonstrate later through 
our finite size spectrum analysis, by duality, the S-matrix of these kinks is 
identical to that of the local spin flip excitations on the paramagnetic side at a 
corresponding coupling, g — > l/g > 1. 

On a ring, PBC implies that 0\ + 9 2 = 0. Furthermore, in contrast to the 
paramagnetic phase, in the ferromagnetic phase one must take into account 
the presence of the three possible vacuum states when constructing periodic 
solutions. A way to do that is by keeping track of the vacuum polarization 
at position x = 0, e.g. As a consequence, wave function amplitudes must also 



have a vacuum label on the ring, As 



A 



</<) 



However, there is a subtle 



difference between scattering in an infinite system and scattering on the ring. 
As illustrated in fig. 13.11 moving one of the kinks around results not only in 
a phase change and a collision of the elementary excitations, but the domain 
orientations also change in a peculiar manner: the configuration essentially turns 
"inside out" . Correspondingly, the amplitudes of the wave functions change as 



.4 



e ifcli si(fc! 



hi) A 



(m-i) 



s 2 {k 1 - k 2 ) A { 1 



+ 



(17) 



The states discussed so far are not eigenstates of the cyclic operator, Z (cf. 
eq. ©). However, we can define eigenstates of Z by taking linear combinations 
of them. Combining e.g. the three ferromagnetic ground states we find 



IQ) 



1 

71 



-iQQfj, 



10), 



(18) 



Similarly, we can define the two particle states, \x\ 6%, x 2 6%; Q), and the cor- 
responding scattering states and wave function amplitudes, A^ e ^ by simply 
mixing the states and the wave function amplitudes as in eq. (|18[) . Since the 
quantum number Q is conserved, the periodicity condition of two-particle states 
simplifies in this basis. Taking the kink i = 1 around the ring, the relation 
in eq. ([T7|) implies the following equation for the amplitudes C\ = A ( +_ and 

C 2 EE A^ 



-+ ' 



Ci 
C 2 



Ak\ L 



si(ki - k 2 )e- in ® s 2 (h ~ k 2 )e~ in ® 



s 2 (h - k 



2)e 



s\{ki 



C 2 



(19) 



and a similar equation is obtained for moving around particle i = 2. The 
structure of these equations is analogous to those in the paramagnetic case, but 
the scattering lengths are replaced by some effective Q-dependent scattering 
lengths, a t , s —t bf s , yielding the triplet and singlet spectra 



E%% 2 -2A = £ 



0(l/L 2 ) 



1 2 1 (m - n 2 + l) 2 

-(ni+n 2 ) +- -g 

4 4 (1 + 1^)2 



4>i+-) 2 + 4^ + ^ 2 ) 



(20) 



Here the lengths bf and bf can be obtained by expanding the phases of the 
eigenvalues of the matrix in eq. (fHfl) , 



s% t {k) ee e 2i O fe ) = Sl (k)cosnQT Y / sl(A;)-sf(/c)sin 2 VlQ , (21) 
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for low momenta. In the Q = sector the scattering lengths are thus given by 
bf~° = a t and b® =0 = a s . The finite size spectrum is thus in agreement with 
the duality relation, eq. (J3J), provided that 

a t {g) = a t (l/g) , a s {g) = 0.(1/5) . (22) 

In the Q = ±1 sector we get, on the other hand, 

6? =±1 = ^-^,&? =±1 = ^-^f (23) 

Equations (|20|) and (|23|) arc the most important predictions of the effective 
theory. Together with the duality relation between the ferromagnetic and para- 
magnetic phases, they allow us to fit the DMRG data on the ferromagnetic side 
g < 1 without any free fitting parameter. 

4. Numerical results 

4--1- Technical details 

In the numerical calculations, we find the lowest lying eigenstates of the 
Hamiltonian eq. ([1} using the lattice units J = a = H = 1. We perform a 
standard DMRG calculation [25], where we make use of the §3 symmetry in the 
Q sector to perform a sub-blocking of the vector space. We use a two site A»»B 
super block configuration and targeted for up to 15 states lowest in energy. In 
order to achieve convergence for the larger system sizes, wc start with an initial 
run of targeting the lowest 3 states only, performing 11 finite sweeps and keeping 
2000 states per A/B block. Wc then restart this run increasing the number of 
low lying target states to 7 and continue with 11, 15 low lying states keeping 
2500, 3000, and 4000 states per A/B block performing 5 finite lattice sweeps 
in each restart. For the 240 site systems we continued restarting the DMRG 
runs keeping the 15 states lowest in energy and up to 5000 . . . 7000 states per 
block. In order to deal with the degeneracies and the large number of low lying 
states we use the generalized Davidson algorithm ensuring that the solution of 
the preconditioncr u n is orthogonal to the previously found state it„_i which 
helps to avoid stagnation of the Davidson algorithm without paying the full 
overhead of a Jacobi-Davidson scheme. Calculations were performed with a 
multi threaded code running on eight core machines with 64GB of RAM. 

4-2, Single-particle levels 

The numerically implemented PBC forbids single quasiparticlc excitations 
on the ferromagnetic side, where they can appear only under twisted boundary 
conditions^ In contrast, in the paramagnetic phase the charge of quasiparticles 
is Q = ±. Therefore, while we could not investigate single quasiparticle exci- 
tations in the ferromagnetic phase, we could study them in the paramagnetic 
phase in the Q = ± sectors, where they appear as the lowest-lying excitations. 



3 Domain wall excitations of the ferromagnetic phase have a S3 charge Q = 0. 
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Figure 3: (Color online) Single-particle parameters extracted from the DMRG data. Circles 
denote the quasiparticlc gap A(g), while diamonds correspond to the quasiparticle mass m(g). 
Units ofJ = a = h= l are used. The standard error estimated for the fitting is less than 
the linewidth. Inset: quasiparticle dispersion relation e n — A as a function of the momentum 
k n = 2irn/L. 

Equation (Q| and PBC imply in the Q = ± sectors of the paramagnetic phase 
that, for very large systems, the single-particle energies are given by 

e„( ff >l) = A( ff) + ^(fn) 2 + ..., (24) 

with The quasiparticlc gap can thus be identified as 

A(g > 1) = lim (£l = o (.9, L) - E^=° (g, L)) , (25) 

and can be obtained from extrapolating the corresponding numerical data to 
L _1 — > 0. The quasiparticle mass m{g) can be defined and extracted through a 
similar extrapolation procedure. 

The single-particle parameters obtained this way are shown in fig. [21 The 
inset demonstrates that the quadratic dispersion is indeed consistent with the 
numerically computed excitation spectrum. Both m and A decrease as the 
coupling approaches the critical value, g — > 1, where the gap is supposed to 
vanish as A ~ J\g — 1 1 5 ^ 6 - The data are consistent with this power-law behavior, 
but it is difficult to extract the precise value of the critical exponent from them. 

The fact that due to the lack of the single-particle excitations, we cannot 
obtain the quasiparticlc parameters on the ferromagnetic side, g < 1 , directly 
from the DMRG data is of little concern. The duality relation, eq. ([3]) relates 
the quasiparticlc gaps and masses in the two phases, since for two remote quasi- 
particlcs in a very large system we must have 

e(h,g) + e(k 2 ,g)=g (e(ki, 1/g) + e(k 2 , l/g)) . (26) 

This can hold for all momenta k\ , k 2 only if 

m(g) = - m(l/g) , A(g) = g A(l/g) . (27) 
9 
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a) 



£ ni,n 2 


(ni,n 2 ) 


parity 


DMRG 


1/4 


(0,0) 


t 


0.21995 (lx) 


1/2 


(1,0),(0,-1) 


s 


0.50658 (2x) 


1 


(1,-1) 


s 


1.02647 (lx) 


5/4 


(1,0),(0,-1) 


t 


1.13124 (2x) 


5/4 


(1,1),(-1,-1) 


t 


1.21941 (2x) 



b) 



diag 

"1,12 


(ni,n 2 ) 


parity 





(0,0) 


t 


1/2 


(1,0), (0,-1) 


s 


1/2 


(1,0),(0,-1) 


t 


1 


(1,-1) 


s 


1 


(1,-1) 


t 


1 


(1,1),(-1,-1) 


t 



Table 1: a) Asymptotic values of normalized two-particle energies, e ni ,n 2 = (En 1 ,n 2 ~ 
2A)/£o, for L — > oo as predicted by the effective Hamiltonian for a reflective S-matrix in the 
Q = sector of the paramagnetic phase, and the corresponding rescaled energy values, from 
DMRG (without extrapolation to L — > oo). The DMRG data are taken at g = 2 for L = 240. 
b) asymptotic values of energy levels assuming a diagonal S-matrix. 



4-2.1. Two-particle levels: paramagnetic phase 

The prediction of our effective field theory is that for very large system 
sizes, L — > oo, the excitation spectrum becomes universal in the sense that the 
rescaled energies, e ni .n 2 = (E ni , ri2 — 2A)/£rj, approach universal fractions and 
have corresponding universal degeneracies. Furthermore, eqns. (|13[) .([T5 |l and 
(|20| also predict that corrections to this universal spectrum can be fitted just in 
terms of two scattering lengths, a s and at for all levels. The predicted universal 
spectrum, and its comparison with the numerically obtained finite size spectrum 
is shown in Table [TJa. Already without incorporating finite size corrections, 
a very good agreement is found: all degeneracies as well as the approximate 
energies of the states agree very well with the predictions of eqns. ([13)) and (|T5l) . 
However, as we demonstrate in Table [TJb, the finite size spectrum is completely 
inconsistent with the spectrum associated with a diagonal S-matrix: there the 
phase shifts <5^iag an( l ^diag vanish for small momenta, and the asymptotic values 
of the normalized levels are given by (ni + tt^) 2 /4 + (ni — ri2) 2 /4, with n\ > ri2 
and ri\ > n 2 for the triplet and singlet sectors, respectively. We remark that 
the same (inconsistent) values are given by the perturbed CFT calculations, 
discussed in Section 

An even more consistent picture based on the effective theory is obtained 
if one also incorporates 1/L corrections due to the finite scattering lengths, a s 
and at- The latter quantities can be extracted from the finite size spectrum by 
using only the lowest two excited states in the Q = subspace 



L 

at[g > 1) = - hm - 



a s (g > 1) 



lim — 

L— ¥ oo 2 



(E^(g,L)-2A(g))/£ (g,L)-i 
(E^(g,L)-2A(g))/£ (g,L)-\ 



(28) 
(29) 



respectively. The (/-dependence of the extracted scattering lengths, a s (g) and 
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Figure 4: (Color online) Scattering lengths ctt(g) and a s (g) as a function of g in units of the 
lattice constant, a. Inset: the logarithms of the scattering lengths as functions of ln(g — 1) 
show a power law divergence. 



a t (g) is shown in fig. 3J Similar to the correlation length, a s and a t both seem 
to diverge at the critical point, but it is not possible to extract an accurate 
exponent from our numerical data. 

Although the extrapolations in the case of the scattering lengths are less 
accurate than in the case of the single-particle parameters, we can now plot the 
higher energy levels in the Q = sector using eqns. (fT3)l and (fT5|) . and compare 
them to the appropriately rescaled DMRG data for various system sizes. We 
find a very good agreement between the numerical results and the predictions for 
the spectrum of the effective model, as can be seen in fig. [5) A clear convergence 
to the asymptotic values is observed for large L's, and deviations appear only at 
smaller system sizes or at higher energy levels, where the asymptotic description 
must break down. 

4-2.2. Two-particle levels: ferromagnetic phase 

As we discussed before, the spectrum on the ferromagnetic side g < 1 with 
PBC does not have any single-particle levels, from which we could get the quasi- 
particle parameters directly. Although, in principle, it would be possible to fit 
the quasiparticle gap and quasiparticle mass along with the scattering lengths 
from the two-body spectra given by eqns. (f2"U)) . this is not needed. As discussed 
earlier, the duality relation, eq. connects energy scales for g <-> 1/g, and 
thus A(<?) and m(g) through eq. ((27)) in both phases. In addition, it also im- 
plies that all length scales emerging in the problem must be invariant under the 
duality transformation jf) 1/g: relevant length scales appear in the finite size 
spectrum as cross-over scales, and by the invariance of the spectrum, they must 
transform similar to the scattering lengths, eq. (|2"2"|) . 

Fig. [5] provides an explicit numerical evidence for the duality relation in 
the Q = sector. There we show that the appropriately rescaled DMRG data 
obtained for g < 1 completely overlap with the data on the paramagnetic side 
g — > 1/g. This gives a numerical proof for the relations ([3]), (f2"2"j). and (|27p . We 
emphasize that the duality relation holds for all system sizes L. 

In addition to the consistency between the effective theory with a reflective 
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Figure 5: (Color online) Comparison of the rescaled energies from DMRG and the effective 
theory for g = 3, g = 2, and g = 3/2 in the Q = sector. The scattering lengths at and 
a B were fitted using the lowest lying two levels, respectively. No further fitting for the higher 
levels was used. Using the duality of the model, we also show the rescaled ferromagnetic 
spectrum for g = 1/g in the Q = sector. 
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asymptotic S-matrix and the numerical data in the Q = sector, probably the 
most important check of validity of our assumptions is the comparison to the 
DMRG results in the ferromagnetic Q = 1 sector. There the predicted finite size 
spectrum is given by eqns. ([20)1 together with eq. ([23"]) . We emphasize that all 
parameters of the effective theory are fixed already and no further adjustment 
to the theoretical spectra is possible. As we show in fig. [51 the effective theory 
presented here also describes the numerical data in the Q = 1 sector within 
numerical precision in the regime, £ <C L, and fully confirms eqns. (|20[) and 

(ESI. 



5. Comparison to the scaling Potts field theory 

5.1. Scaling Potts field theory as a perturbed conformal field theory 

Here we only give a brief review of the scaling Potts field theory; a detailed 
analysis of the scattering theory is given in a separate paper [l8| . The scaling 
limit of Potts model at the critical point is a minimal conformal field theory 
with central charge C = | 0, Q . The Kac table of conformal weights is 



{h r , s } 



( I ! f 3 \ 

| f8 h 5 ¥ * 

V a T 3 8 U / 



r = l,...,4; 

8 = 1,..., 5. 



The sectors of the Hilbert space are products of the irreducible representations 
of the left and right moving Virasoro algebras which can be specified by giving 
their left and right conformal weights as 

S h ,h = Vh 8> V- h . 

There are two possible conformal field theory partition functions for this value 
of the central charge Q ■ The one describing the Potts model is the D4 modular 
invariant, for which the complete Hilbert space is 

u = s , e S2 2 © S? i © s 3 3 

5 1 5 5 ' 5 

_l_ S 1 1 S~2 2 $2 2 

15 ' 15 15 ' 15 3 ' 3 3 ' 3 

®Siz ©5^2 ©S ,3 ©S 3 ,o • ( 30 ) 

Note that not all of the possible representations occur in the Hilbert space; 
there is another modular invariant partition function called A4 which includes 
all sectors of diagonal form Sh,h allowed by the Kac table exactly once: % = 
©I, V/j <X> Vh ■ The A4 model corresponds to the scaling limit of a higher multi- 
critical Ising class fixed point with symmetry Z2. In contrast, the D4 conformal 
field theory is invariant under the permutation group S3 generated by two ele- 
ments Z and C with the relations 

Z 3 = 1 , C 2 = 1 , CZC = Z^ 1 , 

which have the signatures sign Z = +1 and sign C = — 1. The sectors in the 
first line of (|30[) are invariant under the action of the permutation group, S3, 
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Figure 6: (Color online) The spectra for g = 1/3, g = 1/2, g = 2/3 in the Q = 1 sector. Only 
the numerically converged data points and the corresponding curves from the effective theory 
are shown. 
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while the two pairs on the second line each form two-dimensional irreducible 
representations, as characterized by the following action of the generators: 

C|±)HT>, Z|±)=e ± ^|±). (31) 

Finally, sectors in the third line of cq. (|30[) form onc-dimcnsional signature repre- 
sentations, where each element is represented by its signature. These sectors are 
in one-to-one correspondence with the families of conformal fields: the primary 
field in the family corresponding to S h } L has left and right conformal weights h 
and h, and a corresponding scaling dimension, A h j L = h + h, and is denoted by 
<5> h ^, with an optional upper ± index for fields forming a doublet of §3. 

Relevant fields are exactly those for which h + h < 2. It is then obvious 
that the only Z3-invariant spinlcss relevant field is $22, which means that 
the Hamiltonian of the scaling limit of the off-critical Potts model is uniquely 
determined 6] . 



H = H Sf +r J dx $22 . (32) 

The sign of the coupling constant r corresponds to the two phases: r > is 
the paramagnetic, while r < is the ferromagnetic phase. Up to normalization 
factor, it is given by 

r <x (g - l)a- 6 / 5 , (33) 

with a the lattice spacing. The scaling limit is achieved by taking a — > and 
g — > g c = 1 such that r remains finite. In this limit, the gap 

A^ r 5/6 „ | 5 _l| 5 /6 hc/a 

remains also finite. Here, for clarity, we restored H and c, which are usually both 
set to unity in relativistic quantum field theory. 

The scaling Potts field theory (|3"2")l is known to be integrable (l3| . and its 
spectrum and scattering matrix was determined exactly (l3 . Il6l ] . In the para- 
magnetic phase, the vacuum is non-degenerate and the spectrum consists of a 



pair of particles A and A of mass m, which form a doublet under Z3 [151 ]: 

C\A(P)) = \A(P)), Z\A((3))=e iQ \A((3)), 

C\A(P)) = \A((3)) , Z\A(f3)) = e- iQ \A(l3)) . (34) 

The excitations A and A correspond to the local spin flip excitations of chirality 
a = Q = ± of the lattice. The generator C is identical to charge conjugation (^4 
is the antiparticle of A). Choosing units in which h — c = 1, two-dimensional 
Lorentz invariance implies that the energy and momentum of the particles can 
be parameterized by the rapidity j3: 

E = m cosh [3, p = m sinh ft . 



The two-particle scattering amplitudes are 13 1 



sinh (8%. + f) 

s AA {fti2) = s AA (ft 12 ) = — y — 4 

sinh {§12. + ZLi N 

S AA (fti2) = S AA (ft 12 ) = ^ ^f, (35) 

sinh ( £f - r 
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where /?i2 = /?i — /?2 is the rapidity difference of the incoming particles. This 



S matrix was confirmed by thermodynamic Bethe Ansatz [261 ]. We remark 



that the pole in the Saa = Saa amplitudes at /3i 2 = ^ corresponds to the 
interpretation of the particle A as a bound state of two particles .4, and similarly, 
A as a bound state of two ^4s, under the bootstrap principle (a.k.a. "nuclear 
democracy"). The pole in S a a = $AA amplitudes at /3i2 = ^ has a similar 
intcrpretation in the crossed channel, and it does not correspond to a true bound 
state in the neutral sector. 

The excitations in the ferromagnetic phase are topologically charged [l6j |. 
Similar to the lattice model, the vacuum is three- fold degenerate |0) M (/i = 
— 1, 0, 1). The action of S3 on the vacua is 

£|0)„ = |0Wlmod3, C|0) M = |0)_ M , 

and the excitations are kinks of mass m interpolating between adjacent vacua, 
and correspond to domain walls on the lattice. The kink of rapidity f3, interpo- 
lating from /i to // is denoted by 

K^i (/?) , fj, - /jf = ±1 mod 3 . 

The scattering processes of the kinks are of the form 

K^ifc) + K Vil ,(fa) K^ifii) + K vltl ,(fo) , 

with the scattering amplitudes equal to 



S [ V v V ] CS12) = \ n ro , , r , (36) 



{Saa(Pi2) if M =m' 



This essentially means that, apart from the restriction of kink succession dic- 
tated by the vacuum indices (adjacency rules), the following identifications can 
be made 

K- (R\-I A ^ M-^ = +lmod3, (w 
K ^ )= \A(P) „-„ = -! mod 3. (3?) 

in all other relevant physical aspects (such as e.g. the bound state interpretation 
given above). 

The validity of the S-matrix expressions (|3"5T) and ([3"B1 for the scaling Potts 
model can be checked by comparing the finite size spectrum of the corresponding 
Bethe Ansatz equations to that of ((32)) as obtained by the truncated confornial 
space approach (TCSAjl This is performed in detail in [18|. In the TCSA, 



one determines the finite size spectrum of eq. (|32p numerically by truncating 
the finite volume Hilbert space by imposing an upper cutoff in the eigenvalue 
of the conformal Hamiltonian. For the ground state, this is equivalent to the 
standard variational calculus in quantum theory, where the variational wave 
function Ansatz is expressed as a linear combination of a finite subset of the 
cigenstates of the conformal Hamiltonian. By looking at the conformal fusion 



The truncated conformal space approach was originally developed in [27 
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rules implied by the three-point couplings [H, |29|, Hoj], it turns out that the 
perturbing operator acts separately in the following four sectors: 



H = S ,o ® Sl 1 ® Sl 1 ® 5 3,3 , 

5 ' 5 5 ' 5 

n± = si t © st 2 , 

15 ' 15 3 ' 3 

Hi = 5| j ©57 2 ©5 ,3©5 3 ,o , (38) 

so the Hamiltonian can be diagonalized separately in each of them. In the 
lattice language, H± correspond to the sectors Q = ±, while Ho © Hi span the 
Q = sector. Charge conjugation C implies that the Hamiltonian is exactly 
identical in the sectors 1-L+ and Furthermore, the spectrum is invariant 

under transformation, r — > —t in sectors Ho and Hi- This is the consequence 
of a Z2 symmetry in these sectors, which leaves the fixed point Hamiltonian 
and the conformal fusion rules in these sectors invariant^. Away from the critical 
point, it can be interpreted as the continuum form of the duality transformation 
Q in the scaling limit. 



As further discussed in [18j, the detailed TCSA calculations indeed confirm 
that the S'-matrices (|55|) and (ptfj|) correctly describe the paramagnetic and fer- 
romagnetic phases of the scaling field theory. However, as we discuss below, 
the Bethe Ansatz spectra computed with (|35[) and (|36|) both turn out to be 
inconsistent with the numerically computed finite size spectra. 

5.2. Comparing the scaling field theory to DMRG 

In order to compare the DMRG to the scattering matrices ()35I36|) directly, 
we need to rescale the variables to appropriate units in which c = 1. The 
rclativistic relation 

A = mc 2 



allows to determine the speed of light c = y/ A/m in lattice units (aJ/ti). We 
recall that, according to eq. @, A is the infinite volume limit of the energy 
gap between the stationary one-particle state and the ground state, while m can 
be determined from the large volume behavior of the first excited one-particle 
state. We then introduce the dimcnsionless volume variable (h = 1) 

/ = mcL , 

i.e. we measure the volume in units of the Compton length. After rescaling the 
DMRG spectrum to these units, we expect the spectrum of one-particle states 
to follow the relativistic dispersion, 



A (E(L)-E (L)) = yi+^j +0(e 
p/mc — 2it n/l , 

where Eq(L) denotes the ground state energy up to exponential finite size correc- 
tions. [f| The dispersion above indeed describes the numerically obtained finite 



5 The conformal fusion rules do not allow the extension of this symmetry to the H±. 
6 These corrections are due to vacuum polarization and particle self-energy corrections 
induced by finite volume [iilll . 
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size spectrum of low energy quasiparticles, as demonstrated in fig. [7] High en- 
ergy deviations are mainly cut-off effects due to the fact that the DMRG data 
are not close enough to the fixed point. 

Scaling field theory predicts that (neutral) two-particle states in the param- 
agnetic phase are described by the Bethe-Yang quantization conditions 

e a sinh ft _ ^ ) = 1 j Q U sinh P*S AA (P2-Pl) = 1 , 

or, in logarithmic form, 

I sinh/3i + 2S AA {/3 1 - /3 2 ) = 27rni , 
I sinh f3 2 + 28 AA (p 2 ~ Pi ) = 2nn 2 , 

with n\ and n 2 integer quantum numbers, and the phase-shift function defined 

as 

Saa{P) = -\^S aa {P) . 

The Bethe-Yang equations are nothing else than the conditions (fTTj) stated in 
terms of the notations of the scaling field theory. The energy relative to the 
ground state can be computed as 



E(L) - E (L) = A(cosh/3i + cosh/3 2 ) + O (e 



and is accurate to all orders in l/l, as for the one-particle states. To conform 
with the conventions used previously, we plot the rescaled quantity 

P E(L) -2A- E Q (L) 



{2-Kf 



against I, The results are shown in fig. [HI which shows that the scaling field 
theory correctly describes the singlet levels. In the language of the field theory, 
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Figure 8: (Color online) Comparison of rescaled DMRG data (for g = 3) to the two-particle 
levels predicted by the bootstrap S matrix. While the singlet spectrum is perfectly reproduced 
without further fitting parameter, the triplet sector cannot be fitted. 



these are exactly the two-particle levels in the C-odd sector Hi. However, the 
triplet levels cannot be explained by the bootstrap S matrix. In the scaling held 
theory, these levels are in the C-even sector "Ho and are described by the same 
Bethe- Yang equations, which means that they are exponentially degenerate with 
their sin glet counter-parts. While this picture is fully confirmed by the TCSA 
analysis [18[ , it is clearly not consistent with the DMRG spectrum. 

To see the problem more clearly, one can perform a direct comparison of 
the phase-shift function to the DMRG spectrum. Provided the energy levels 
E(L) are known, the Bethe- Yang equations can be used to extract the phase- 
shift function from them. The results are shown in fig. [9] For the singlet levels 
the slope of the phase-shift around the origin agrees quite well with DMRG 
data, which means that the bootstrap S matrix gives correctly not only the 
low-energy value of the phase-shift, but also the scattering length. For larger 
j3 the deviations are explained by cut-off effects since these correspond to lower 
values of the volume, closer to the scale of the lattice spacing. 

However, the phase-shift extracted from the triplet states does not agree with 
the bootstrap prediction at all: neither the low-energy value nor the scattering 
length is consistent as figure [S](b) demonstrates. 



6. Conclusions 

In this work, we presented an effective low-energy theory to describe the 
lattice version of the q = 3 state quantum Potts model in one dimension. In 
our approach, the interaction of the quasiparticles is characterized by a scat- 
tering matrix, which takes the universal exchange form S(k — > 0) — > —X in 
the asymptotic low-momentum limit. A scattering matrix of this form has been 
first motivated in Ref. Q by perturbative and renormalization group arguments. 
To prove (or disprove) the validity of these statements, we calculated the finite 
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Dashed lines are the theoretical predictions from the bootstrap. 
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size spectra corresponding to this effective theory and compared them to the 
" numerically exact" energies obtained by density matrix renormalization group 
calculations as well as to the results from the scaling Potts field theory. The 
results of this comparison are very clear: the numerically obtained spectrum 
fully supports the picture of Ref. |3j , and it rules out the asymptotic diagonal 
S-matrix predicted by the scaling Potts field theory We also find that the boot- 
strap S matrix is able to reproduce only part of the finite size spectrum (singlet 
part), and is thus inconsistent with the numerical data. 

The structure of the S-matrix can be understood on physical grounds. It 
is trivial to show that for interacting bosonic quasiparticles in the absence of 
internal quantum numbers any local interaction is relevant fioj , and leads to a 
scattering phase shift, S = ±7r/2. However, for particles with some nontrivial 
internal quantum number the situation is somewhat more subtle. In the anti- 
symmetric (singlet) channel the wave function has a node, and therefore local 
interactions are asymptotically irrelevant, leading to S s (k = 0) = 0. In contrast, 
in the symmetric (triplet) channel the orbital wave function does not vanish 
when the two particles approach each other. Here interactions are relevant, and 
lead to 6t(k = 0) = ±7r/2. These arguments immediately imply an S-matrix of 
the exchange form. Note that in order to obtain a diagonal S-matrix one needs 
5 t (k = 0) = 0. However, this is only possible for very special effective interac- 
tions between the quasiparticles 0] , and is not guaranteed by § 3 symmetry. 

Our numerical calculations also allowed us to check the duality between the 
ferromagnetic and the paramagnetic states. We have shown that in the Q = 
sector the whole finite size spectrum obeys duality. As a consequence, the scat- 
tering lengths, the quasiparticle masses and gaps also satisfy duality relations, as 
demonstrated by the analysis of the finite size spectra. The extracted scattering 
lengths seem to diverge at the critical point. 

We should also remark that our calculations do not support the existence 
of bound states in the {H — } and { — h} channels: every state in the finite 
size spectrum could be identified as an extended two-particle state. Bound 
states were also absent in the paramagnetic Q = ± sectors, where we have just 
observed extended single particle excitations at and slightly above the gap, A, 
all in agreement with the simple effective theory. We remark that this also 
agrees with the bootstrap, where particles A and A can be interpreted as bound 
states of AA or AA in the spirit of Chew's "nuclear democracy" ; therefore, one 
does not expect any additional state besides the multi-particle states built from 
A and A (or the corresponding kinks in the ferromagnetic phase). 

Let us finally comment on the discrepancy between the scaling Potts field 
theory and the properties of the Potts model on a lattice. We believe that the 
main source of this discrepancy is related to the assumption of integr ability. Al- 
though the q = 3 lattice Potts model is intcgrable at the critical point, g = 1, it 
is believed to be non-intcgrable for other values of g. However, the scaling Potts 
field theory is intcgrable, and the bootstrap, together with the assumption of S3 
symmetry, leads to a diagonal S-matrix for q = 30 Also, in the perturbed CFT 
description, one perturbs the continuum Potts model with the leading relevant 
operator, which leaves the model integrable, and gives a spectrum consistent 



7 Remarkably, for any other q 7^ 3 the asymptotic S-matrix assumes the universal exchange 
form. 
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with the bootstrap S- matrix (as shown in detail by the analysis in [18|). These 
methods, at least in their original form, are thus unable to describe the lattice 
Potts model for j ^ 1, probably because requiring integrability conflicts with 
the true non-integrable nature of the 1+1-dimensional q = 3 state lattice Potts 
model. We think that to describe the q = 3 lattice quantum Potts model, one 
needs to allow for perturbations or cut-off schemes which violate integrability. 
In perturbed CFT, a possible candidate would be adding the leading irrelevant 
operator which would violate integrability. 

Also, we have to admit that our lattice calculations arc still relatively far 
away from the critical point itself. It is very hard to believe that the asymptotic 
theory discussed here would suddenly break down as one approaches the critical 
point, 3=1. In this respect, we think that our results are conclusive. However, 
due to the vanishing of the gap, and the divergence of the correlation length 
and the scattering lengths, the effective theory will clearly be limited to smaller 
and smaller momenta, and correspondingly, to smaller and smaller temperature 
regimes. Indeed, fig. [9jb seems to support the emergence of a new low-energy 
scale, A* ~ h 2 / a 2 (g)m(g) , which vanishes faster than A as g — > g c , in agreement 
with the scenario of irrelevant operators, discussed before. The asymptotic the- 
ory is restricted to be valid below this scale, i.e., for quasiparticles of momenta 
smaller than ~ l/a t . 

Based on the present data we can neither exclude nor admit the possibility 
that at some intermediate energies a description in terms of the diagonal S- 
matrix may be adequate. While in the singlet sector, shown in fig. |H] (a), the 
deviation from the bootstrap solution for intermediate rapidities /3 ~ 1 can be 
explained by the finite UV cutoff, the situation is not so clear in the triplet case, 
shown in fig. [9] (b). On the one hand, as g approaches the critical point, the 
phase shift for a small but fixed value of f3 ^ seems to slowly approach the 
curve corresponding to the diagonal S-matrix. On the other hand, while the 
phase shift seems to have reached the scaling limit for j3 « 1, the deviation from 
the bootstrap prediction remains considerably larger than in the singlet case. 
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Appendix A. Duality 

We show that a one-to-one correspondence (duality relation) exits between 
the energies in the Q = subspace with PBC for couplings g -O- 1/g. To do 
this, we introduce two sets of basis states. The first set is defined using the local 
"spin-flip" states 

L 

\{Xi})=H\^)i, (A.l) 

i=l 

where each A; <G {0, 1, —1}. Restriction to the Q = subspace implies Y^,i=i ^» = 
0. For a given sequence of {A^}, there exists another orthogonal set of dual 
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states, defined by using the Xi as domain wall labels, 

L 



KM) -H{Ai}>M = I 



IM=H + *^2X- \ . (A.2) 

3 = 1 I A 



By construction, (since Yli=\ Ai = 0), these states automatically satisfy PBC, 
but they are not eigenstates of the permutation operator C. However, we can 
construct states within the Q = subspace by defining 

\{Xi})^^[l + Z + Z 2 }\{Xi})^. (A.3) 

Straightforward algebraic manipulation yields that the matrix elements of 
the two terms of the Hamiltonian, Hi = J^i ^t^t+i ano ^ ^2 = J2i Pi m cc l- © 
satisfy the following identities: 

Let us now assume that the states 

I ?1 > = E^{mI{M> ( a - 5 ) 

are normalized, orthogonal eigenstates of H± + gH 2 in the Q = subspace, 

(m.\(Hi+gH 2 )\n) = S nm E n {g). (A.6) 
Then let us define the set of dual states as 

W=E A {A,}KM>- (A.7) 

Using eq. (|A.4[) , we immediately see that that these diagonalize the dual Hamil- 
tonian, Hi + l/gi?2, 

H(#i + 1/.9#2)R) = (l/<?)(n|.gff 2 + Hi\m) = (1/ g)S nm E n (g) , (A.8) 
yielding the duality relation, eq. (|3]). 
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